Testing for Multipartite Quantum Nonlocality Using Functional Bell Inequalities 
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We show that arbitrary functions of continuous variables, e.g. position and momentum, 
can be used to generate tests that distinguish quantum theory from local hidden variable 
theories. By optimising these functions, we obtain more robust violations of local causality 
than obtained previously. We analytically calculate the optimal function and include the 
effect of nonideal detectors and noise, revealing that optimized functional inequalities are 
resistant to standard forms of decoherence. These inequalities could allow a loophole-free 
Bell test with efficient homodyne detection. 



Bell famously showed that the predictions of quantum 
mechanics (QM) are not always compatible with local 
hidden variable theories (LHV) jlj. Surprisingly, this 
fundamental result, which underpins the field of quantum 
information, has not been rigorously tested 0|. There are 
no experiments yet that can eliminate all LHV, either due 
to low detection efficiencies [3, 0] or lack of causal sepa- 
ration. Rigorous tests are also needed to fully implement 
some quantum information protocols, like that of Ekert 
[3] which employs a Bell inequality (BI) as a test of se- 
curity in a cryptographic scheme. All of these early tests 
and protocols employed quantum measurements with dis- 
crete outcomes of spin or particle number. 

In this Letter, we develop functional moment inequal- 
ities to test for quantum nonlocality. We then use varia- 
tional calculus to optimize the choice of measured func- 
tion. As a result, we obtain Bell nonlocality for larger 
losses and for greater degrees of decoherence than pos- 
sible previously. The outcome can be feasibly tested in 
the laboratory, since the detectors required are efficient 
quadrature detectors. More generally, functional nonlo- 
cality measures could lead to new applications in quan- 
tum information. The important advantage is a much 
greater robustness to noise and loss. 

As well as potentially overcoming the loophole problem 
mentioned above, formalisms to test LHV for continuous 
variables provide an opportunity for testing QM in new 
environments, and give a better understanding of the ori- 
gin of the nonlocal features of QM. This is particularly 
true given that entanglement alone does not guarantee 
failure of LHV for mixed states Q . 

With this objective, there is the fundamental question 
of how to quantify the strength of nonlocality, in the ab- 
sence of a single test for nonlocality that is necessary 
and sufficient for any quantum state. Mermin [§] used as 
a measure the deviation of the QM prediction from the 
LHV bound, based on a particular BI. A second strat- 
egy discussed recently by Cabello et al. Q is to quantify 
the strength of nonlocality by the robustness of the vio- 
lation with respect to a decoherence parameter. In this 
approach one determines the critical efficiency r\ or the 
critical degree of purity p required for a violation. Here, 



we evaluate all three measures to show strong correlations 
between them. 

Recently, Cavalcanti et al. (CFRD) showed 0] that 
Bell inequalities can be derived for the case of observables 
with continuous and unbounded outcomes, like position 
and momentum. This approach is significant in estab- 
lishing that quantum nonlocality does not rely on the 
discreteness of the measurement outcomes. Continuous 
variable (cv) inequalities also provide an avenue to un- 
derstanding how manifestations of quantum nonlocality 
can be manipulated by choice of observable. 

The original CFRD inequality is \(Uk=i( x k + 
ipk))\ 2 < (ILk-iixj. +pl)), where x k , p k are the out- 
comes of two arbitrary measurements, represented in 
QM by observables X k , Pk, at site k [ijj. Where X 
and P are quadrature measurements with canonical po- 
sition and momentum commutation relations, CFRD 
showed that the symmetric state [\q)®n /2^®N /2 + 
|l}«W2|o)®tf/2}y\/2 violates the inequality for N > 10. 
In this case, the states |0), |1) are eigenstates of a) a 
where a = X + iP, so the prediction could in principle be 
tested with photonic Greenberger-Horne-Zeilinger (GHZ) 



states produced in the laboratory [11]]. Note that in the 
above state there are N field modes but only N/2 pho- 
tons. It can be prepared from a A r /2-photon GHZ state 
{\H)® N / 2 + \V)® N / 2 }/V2, where \H), \V) represent hor- 
izontally or vertically polarized single-photon states, by 
passing each photon through a polarizing beam splitter. 
These violations are robust with loss. The critical effi- 
ciency rjcrit required for violation tends to r\ cr i t — > 0.81, 
as N — > 00. Quadrature measurements with local os- 
cillators are highly efficient, with reported efficiencies of 
99%. However, generation losses from mode-matching 
can degrade the experimental efficiency, so 81% is still a 
challenging practical benchmark. 

Instead, we introduce a functional moment Bell in- 
equality by considering arbitrary functions of the out- 
comes at each site. This new approach to nonlocality 
utilizes a general functional optimization of continuous 
variable observables. We find the optimal function that 
maximizes a violation of the inequality for a given effi- 
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ciency r\ and state purity p. We show that the optimal 
function has the form x/(l + enx 2 ), where En is a pa- 
rameter related to N and r\. This gives an inequality 
which is violated by the GHZ states of ((6j) for N > 5. 
The violation increases exponentially with N, while r\ cr i t 
decreases asymptotically to 0.69 for a pure state (with 
p = 1), thus dramatically reducing both the number of 
modes required, and the required efficiency. 

When the functions correspond to a simple binning 
of a cv observable to give binary outcomes [l2|, our in- 
equalities reduce to those of Mermin ]8L We extend the 
analysis of Mermin and Acin et al. [13j, and calculate 
results for homodyne detection for more feasible types of 
state. We find that (rjp 2 ) crlt = 2 (1 ~ 2N '>/ N 7r, which gives a 
critical efficiency for a pure state at large N of f] = 0.79. 

Functional Moment Inequalities. We present a proof of 
the functional moment inequality taking explicit account 
of functions of measurements that can be made at each of 
N spatially separated sites. We denote the measurement 
made on the system at the k-th site by X k , and the 
outcome of the measurement by x k , where 9 represents 
a choice of measurement parameter. Bell's assumption 
that LHV can describe the outcomes implies that the 



measurable moments 



N i 



can be expressed in 



terms of a set of hidden variables A as 

{x°4...x%) = J dXP(X)(x e 1 ) x (xt) x ...(x%) x , (1) 

where (x e k }\ is the average of x e k given a LHV state 
A. Next we construct, for each site k, real functions 
of two observables fk(x k ), g k {x k ), and define the com- 
plex function: Fk — fk(x k ) + igk(x k ). The complex 
moment (FiFz—Fn) can be expressed in terms of real- 
valued expressions of the type (fi(xf) giix^ ).../jv(x^)), 
etc. Of course, fk{x 9 k ) is an observable obtained from 
x k by local post-measurement processing. Eq. (Q]) 
must therefore also be valid for (/(^i)/( x 2) ■ ■ ■ f( x N)) = 
j x dXP(X)(f(xD) x (f(xl)) x ...(f(x%)) x . For an LHV, 
the expectation value of products of the Fk must satisfy: 



(F 1 . . . F N ) = / d\P(X)(F 1 ) x . . . (F N ) X 



(2) 



where (F k ) x = (f k (x e k )) x + i{g k {x 9 k )) x . From ©, the 
following inequality must therefore hold: 

\{F 1 F 2 ...F N )\ 2 < f dXP(X)\{F 1 ) x \ 2 ...\{F N ) x \ 2 . (3) 

Now for any particular value of A, the statistics pre- 
dicted for fk{xk) must have a non-negative variance, i.e., 
(fk(xk)) 2 x < (fk(xk) 2 ) X - Writing JSJ explicitly in terms 
of the fk 's and using this variance inequality we arrive at 
the CFRD inequality with functional moments: 
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We will measure the violation of this inequality by the 
ratio of the left- (LHS) and right-hand sides (RHS). 
Defining the Bell observable B = LHS/RHS, failure 
of LHV is demonstrated when B > 1. In order to get 
stronger violation of local causality, we optimize the func- 
tion of observables by considering 
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(5) 



Here, we consider the class of entangled states 

W = (|0)® r |i)® (A, ~ r) + |l)® r |0)® (JV ~ r) )/V2 . (6) 

Thus r — N corresponds to extreme photon-number- 
correlated states, a superposition of a state with pho- 
tons at all sites and a state with 1 photon at each site. 
Next, we consider how to optimize the function f k and 
gk to generate a robustly violated inequality, including 
losses and noise. 

We use variational calculus to find the optimal func- 
tion using the condition of Eq. ([5]). For simplicity, we 
assume the functions fk and g k are odd. The LHS can 
be maximized by choosing orthogonal angles, while the 
RHS is invariant with angles. We find that 
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where I k = 2je 
(f k ) 2 )]dx, and 7° = fe- 2 * 2 l(f+? + (f k ) 2 )]dx are dif- 
ferent integrals for x which contribute to the expectation 
values in both sides of inequality |(4]). Here f k = f k ±g k , 

9 2 

and the factor e~ lx was obtained from the joint prob- 
ability of observables. Requiring SBN/Sf k = 0, we find 
the optimal condition: fk(x) = ±g k (x). The components 
of complex functions f k , gk are the same at each site, and 
have the form 



fk{x) = g k {x) 



1 + SnX 2 



(8) 



For the even N case, it is optimal to choose r — N/2. 
Then en is independent of TV, but has to be calculated 
numerically since it satisfies a nonlinear integral equa- 
tion: e n = 4J°/I. 

For N an odd number, the greatest violations occur 
for r = (N — l)/2. The optimal function has the same 
form as in J5]) except that the parameter en changes to 



s N , where: 



£ N 



en 



Ne+ 



-N 



Ne+ 



-N 



(9) 



and = En ± 4. However, the numerical value of En 
and e' n now depend on N, as the integral equation ([9]) 
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for odd values of TV is TV-dependent. This provides bet- 
ter violations of than any other arbitrary function, 
provided TV > 5. The maximum B N value with this op- 
timal choice is shown in Fig. HJ compared with the CFRD 
result which uses a simple correlation function. 




Figure 1: Maximum violations of functional cv inequality with 
GHZ states as a function of the number of modes. The vio- 
lations using the optimal function (solid) are much stronger 
than the CFRD result (dashed). 

Binned cv outcomes for Mermin-Klyshko inequality 
(MK) We will also briefly consider binning meth- 

ods. Specifically, we define the binning functions f k (x) = 
9k{x) = fbin{x) = +1 if x > and —1 otherwise. For 
such discrete outcomes, the original formalisms of Mer- 
min and Klyshko @, can be used. The CFRD in- 
equality for discrete outcomes reduces to that of Mermin 
as can be seen by noting that [fun] 2 = 1. Here 
it is known that the Bell inequality introduced for this 
discrete case by Klyshko is stronger. Defining F k = 
fbin(x e k ) + ifbin(x e k ) and Hat = Yl k=1 F k , we can use the 
MK inequality \S N \ < 1 , where S N = 2~ N / 2 [Re{U N } ± 
Im{ITv}], for TV even, and S N = 2-( JV - 1 )/ 2 Rc(Im ){n jv } 
or S N = 2-<- N - 



-l)/2 



V(Re{ITv}) 2 + (Im{n N }) 2 for TV 
odd. These BI have been considered recently for the 
case of extreme photon-number correlated states, where 
r = TV, by Acin et al. We can also define II at by exchang- 
ing the local observables, to obtain similar inequalities. 

We generalize this approach to account for more gen- 
eral angles and states. We find an optimal violation of 
Bn = -^(^) JV / 2 , for arbitrary r and TV, with the optimal 
phases: 6 k = (-l) w+1 7r(fc - 1)/(2N), 6 k = k +n/2 for 
k < r, and 6 k = {-l) N <ir(k - 1)/(2TV), 0' k = 6 k - it/2 
for k > r. This result has been presented by Acin et 
al. for the special case of r = TV. We confirm the 
exponential increase with number of sites TV, but also 
make the observation that the violation occurs for all 
types of states of the form (jfij), independently of r. This 
contrasts with the result for the CRFD inequality, which 
requires r ~ N/2 for violation. However, as explained 
earlier, the states with r = N/2 are straightforwardly 
feasible given a polarisation GHZ state, as opposed to 



the extreme photon-number correlated states considered 
by Acin et al. Therefore, this is a very important exper- 
imental advantage. Violation of the MK inequality with 
binning is possible for N > 3, but, as we will see, this 
strategy is sensitive to losses and noise. 

Sensitivity to loss and state impurity. The value of the 
Bell observable B n increases with the number of sites TV, 
so this is suggestive of a strategy that will allow genuine 
loophole-free violations of local causality. However, it 
may be argued that since increasing the number of sites 
will increase the number of detectors required, there will 
be no advantage. Only careful calculation of the Bell 
observable B n including the detection efficiency rj can 
determine whether the strategy is advantageous. 

Loss is modeled as follows. The field modes a k at each 
site are independently coupled to a second mode a k>vac 
respectively, assumed to be in a vacuum. Photons are 
lost from the field into the vacuum mode, the strength of 
coupling determining the rate of loss. This beam splitter 
model gives the final detected and vacuum mode in terms 
of the inputs a and a vac 



a-out = y/rja + V 1 - V c 

^vac.out 



(10) 



where r\ is the efficiency, the probability of detecting 
a photon after coupling. Since we only measure the 
"a OMt " not the li a va c,out", we need to trace over the latter 
modes to obtain the final density operator p out for the 
detected modes after loss. We can also examine the ef- 
fect of impurity, by considering a state p' = p\ijj)(ip\ + 
(1 — p)p m ix, where p m ix is the mixed state obtained 
with a model for decoherence in the occupation-number 
basis, i.e. Pmix = [|0)® r |l)»( Ar -'-)(0| g)r (l|®( Ar - r ) + 

|l)®r| - > ®(JV-r)^ 1 |®r^ |®(JV-r)]/ 2j and p ig the prob ability 

the system is in the original pure state J6]). 

Including the effect of detection inefficiencies and noise, 
the parameter En is changed to £n(v) f° r the optimum 
function. For TV even we find that 



£at(j?) 



2rj£ 
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B 



N 



2r)+{l-rj)e N ' 
"2(/+) 4 M 2 



,N-2 



irI°C 



(11) 



where £jv is defined as before, and C = r)I+(l— rj)I . For 
the case of odd TV the relevant integral equations change, 
giving a modified (and slightly reduced) Bell variable B' N , 
where: 

Nejf(r]) - e N (T])e N /e N 



(12) 



Here Sj^(rf) = £n(v) + 4, and B n is defined as in Eq (jTTJ) . 
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gled states and application of this method to tests of 
other forms of nonlocality — i.e., entanglement [lH and 
EPR steering [3]. 
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Figure 2: (a) The critical minimum detection efficiency rj cr it 
for pure state, and (b) the critical purity p C rit for ideal detec- 
tors required for violation of functional moment, CFRD and 
MK inequalities with optimal choice of parameters. 



This approach is applied to enable a prediction of 
the effect of loss and noise on the functional inequali- 
ties, and the results are plotted in Fig. [2 These re- 
sults can be compared with the MK binning approach. 
With the choice of optimal angles, we find the values 
of the MK Bell observable with binned cv outcomes is 
B N {p,ri) = ^-(^-) N/2 , which gives the effect of de- 
tection inefficiencies and noise for the optimal choice of 
angles. That implies a critical minimum efficiency and 
purity (rjp 2 ) cr it = 2 < - 1 ^ 2N ^ n tt in order to violate the in- 
equality. For lower N, the strategy of binning and using 
the MK inequality shows an advantage, by allowing a vio- 
lation for N = 3,4, 5, — but even if p = 1, high efficiencies 
?7 > 0.99, 0.93, 0.90 are required. While high detection 
efficiencies are feasible for homodyne detection, these ef- 
ficiency and purity values are still quite challenging once 
generation losses are also taken into account. In view of 
this, the high requirement for r\ cr it for the case N = 3 
may be prohibitive. 

These results show that the functional inequality has 
much greater robustness against noise and inefficiency 
than the MK inequality. For N > 7, the functional cv 
inequality used with an optimal function allows viola- 
tion of LHV at much lower efficiencies and larger max- 
imum noise. The asymptotic decoherence product is 
(pv)oo ~ 0.6918 in the large N limit. For a moderate effi- 
ciency T] cr it ~ 80% one requires N = 10 if the optimized 
function, while the binned MK case requires N ~ 40. 

In conclusion, we have developed a new direction for 
the analysis of cv nonlocality. For the input state treated 
here, the optimal measured function always has the same 
functional form apart from changing the parameter e, but 
more generally, the functional form may depend on the 
experimental decoherence. Future research may include 
further optimization of the functions for different entan- 
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